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Abstract 

An algebraic approach to neutron scattering on a one dimensional potentials is gener- 
alized to diffraction on three dimensional single crystals. 

I hate and will always hate dynamical diffraction 

R.Golub 



1 Introduction 

I feel a solidarity with R.Golub with respect to the standard dynamical diffraction theory, 
because for me it was always so boring to explain diffraction effects by motion of some points 
along mathematical dispersion surfaces. Here I want to present an alternative theory, which 
has a clear physical meaning at every step. My approach is similar to that of |UI2|- 

1.1 Algebra of neutron scattering on two semitransparent mirrors 

To explain the approach, let's consider reflection and transmission of a system of two separated 
semitransparent mirrors. We consider our primary neutron to be described by a plane wave 
tp (k,r,t), and two mirrors — by two different potentials 1*1,2, separated by a distance I as 
shown in fig. ^ If we know reflection 7*1,2 and transmission tip amplitudes of every separate 
mirror, we can find reflection R12 and transmission T12 of both mirrors. 

Suppose the amplitude of the incident wave exp(ikx) on the surface of the first mirror at 
x = is equal to unity at the surface x = of the first mirror. Let's denote X the amplitude 
of the wave incident on the surface of the second mirror (see fig. Q). For X we can construct 
an equation 

X = exp(ikl)ti + exp(ikl)riexp(ikl)r2X. (1) 

The second wave in the right hand side is generated by the X: the incident wave with amplitude 
X reflects from the second mirror, goes to the first one, reflects from the first mirror, and arrives 
to the second one. During propagations between mirrors the wave acquires the phase factor 
exp(ikl). At the second mirror the second term of superposes with the wave, which is 
transmitted through the first mirror and is represented by the first term in the right hand side. 
The sum of two waves according to definition is equal to X. 
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Figure 1: Two semitransparent mirrors are represented by two potentials separated by distance 
I. 



Solution of the Eq. (JTJ) is 
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1 — r\r 2 exp(2ikl) 
Reflected and transmitted waves can now be deduced as 

R12 = r\ + ti exp(ikl)r 2 X, T 12 = t 2 X. 



(2) 



(3) 



Reflected wave is a superposition of two waves. The first one is related to the wave directly 
reflected from the first mirror, and the second one, is related to the wave X, which after 
reflection from the second mirror arrives to the first one and transmits through it. Substitution 
of X into (J2J gives 
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(4) 



Two mirrors above were separated by distance I. It is evident that (@J) are valid for arbitrary 
I. In particular I can be zero, and with such I we have 
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This result is a discovery It shows that an arbitrary potential can at any point be split into 
two parts. If both parts are symmetrical, we get However the potential and its parts, as 
is shown in fig. |2l are in general nonsymmetrical. In that case reflection ri and transmission ti 
amplitudes from the left and from the right are in general different, so formulas JHJ) are to be 
changed as follows 
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where arrow shows the direction of the incident wave 
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Figure 2: Every potential can be split into two parts at an arbitrary point. 



1.2 Properties of reflection and transmission amplitudes 
1.2.1 Nonsymmetric potentials 

To find the difference of amplitudes for waves incident from the left and from the right it is 
sufficient to consider the simplest nonsymmetric potential, which consists of two symmetrical 
ones as shown in fig. El 

Transmission and reflection of the total potential are given by (0). From them we can 
immediately see that T 12 = T 2 i, therefore for arbitrary nonsymmetrical potential we have 




Figure 3: A nonsymmetric potential constructed with the help of two symmetrical ones. 



t = t = t. Reflection amplitudes in general are not equal to each other, however in the case 
of negligible losses they differ only by a phase factor, i.e. *f~ = exp(ix)^ , where x is a re& l 
number. Indeed, because of unitarity we have 

\Y\ 2 + \t\ 2 = | r>| 2 + \t\ 2 = 1, (8) 

therefore \*r~\ = |~r*1. 



1.2.2 Relation between phases of r and t for symmetrical potentials 

The amplitudes r and t of a symmetrical potential can be represented as r = |r| exp(i0. r ) and 
t = \t\ exp(i0t). It is easy to show by a gedanken experiment fl] that for real potentials 

<f) r = (jn± tt/2. (9) 

It follows that t = ±i\t\e l< ^ r , combinations r ± t are unit complex numbers, and r 2 —t 2 = e 2l< ^ r . 

For nonsymmetrical potentials we have 1* = \r\ exp(i^ r ), *r~ = \r\ exp(i^x r ), and instead 
of (0) we can obtain (~^ r + < x r )/2 = Xt i 7r /2- 



1.3 An example of application of the above algebra 

One of the simplest is a rectangular potential shown in fig. 0] If we denote the amplitude 
of wave incident upon the right edge of the rectangle from inside the barrier, we can find an 
equation for X 

X = exp(ik'd)t + exp(ik'd)(— r ) exp(ik'd)(—r )X, (10) 



where 
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and k' = \Jk 2 — u is the particle wave number inside the potential. 
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Figure 4: Rectangular potential barrier of height uq and width d 



The eq. (II Uj) is easily solved. With this X we find both reflection R and transmission T 
amplitudes of rectangle. We omit here all the mathematics (look, for instance, in [21 El) and 
put down the final result 
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which will be used later. 

We restrict here ourselves with only this example of application of above algebra, but those, 
who want to see more examples, can find them in references j3] — jH] 



2 One dimensional periodic potential 



Now we consider one dimensional periodic potentials. Our goal is to find reflection R^ and 
transmission T/v amplitudes of an arbitrary periodic potential, containing N periods, if we know 
respective amplitudes r and t for a single period. However we start with semiinfinite periodic 
potential shown in fig. In fig. |5]we see an incident wave from the left, and the reflected wave 
with an amplitude R. We also see the amplitude X of the wave incident on the second period. 
For X and R we can immediately write the equations 



X = t + rRX, 



R = r + tRX. 



(13) 



For simplicity we assume that the period is symmetric. In writing of this system we used 
and (jHl) for I = and took into account that reflection from the semiinfinite potential 
without one period is the same as from the whole potential. Our system has a very interesting 
symmetry: the equations transform to each other when we exchange r and t. 

We can resolve the first equation with respect to X and substitute it into the second. As a 
result we obtain the equation 

R = r + t 2 R/{l-rR), (14) 

which can be reduced to the form y 2 — 2py +1 = 0, where y = R and p = (r 2 + 1 — t 2 )/2. 
Solution of it is trivial: y = p — i/p 2 — 1 ; however we prefer to represent it as follows 
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and we immediately obtain 
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We know that a shift by one period inside the periodic potential is accompanied by multiplica- 
tion of the wave function by the Bloch factor exp(zga). The amplitude X represents the incident 
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Figure 5: Scattering of a plane wave from a semiinfinite periodic potential. 



amplitude after shifting inside the potential by one period. Therefore X is precisely this Bloch 
phase factor: X = exp(iqa). Because of symmetry of the system (fTHj) we can immediately find 
the Bloch phase factor from ([16)1 . by exchange of r and t: 



[1 + tY -r 2 - J(l-t) 2 -r 2 
X = e iqa = \ — V (17) 
;i + t) 2 -r 2 + J (I - i) 2 - r 2 



There are many modifications of these formulas. We present here only one of them. If we 
represent r = |r|exp(zx r ), and t = ±i\t\ exp(ix r ) , take into account that r 2 — t 2 = exp(2ix r ), 
and substitute into (fTtjj) and (fT7|) . we obtain 



p /CQSXr + \r\ - y/COSXr ~ [j ^ _ ^SUl Xr =F fI ~ V sin Xr ± \t\ 



cosxr + \r\ + Jcosxr - \r\ yshiXr =F \t\ + ysinxr ± \t\ 

The obtained formulas contain full the information about Bragg refections. All these refection 
are total, and the refection amplitude is a unit complex number R = exp(z<f>), when cos 2 \r ~~ 
\r\ 2 < 0. Positions of the Bragg peaks are determined by the equation Xr — (n+ l/2)n, the 
Bragg peaks widths are determined from the equation |A%| < \r\. At the Bragg peaks the 
Bloch wave vector becomes a complex number qa = im + iq"a, which means that the wave 
exponentially decays inside the medium. 



2.1 The finite number of periods N 

To find Rn and Tjv it is sufficient to take into account that a semiinfinite potential with period 
a is also periodic with period Na, and the Bloch phase factor of such a periodic potential is 
X N = exp(iqNa). With this in mind we rewrite the system ()13|) in the form 

X N = Tjv + RnRX n , R = Rn + Tj^RX N , (19) 

where we replaced r and t by i?/v and T/y. The system (fTH|) can be resolved with respect to Rn 
and T/v- As a result we obtain 

^ = i? l g2 , 9 J v T N = exp(iqNa)- — - — . 20 

1 — R z exp(2iqNa) 1 — R z exp(2iqNa) 



These formulas are similar to those for rectangular potential l|12l) . 



2.2 Practical application 

The method demonstrated here can be used for preparation of mirrors with desirable properties. 
In particular, to increase critical angle we need to prepare a supermirror, which consists of 
several periodic chains with different periods and different number of periods jU] . One can also 
prepare a monochromatic mirror, which reflects particles of given energy. For that we need 
only to prepare a periodic system with some special period, giving Bragg diffraction at some 
energy. To exclude Bragg reflection of higher order the single period of the chain must be of a 
special form, which can be well calculated [Hj. 

The method is not restricted to scalar particles. It can be used for neutrons with spin 
in magnetic fields jH], for x-rays, and for elastic waves, and even for investigation of particles 
diffusion and moderation. 

More over it can be used, and we demonstrate it in the next section, for description of 
diffraction in single crystals. With it we have the new dynamical diffraction theory [TO] . 



3 Diffraction in a single crystal 



Interaction of neutrons with media containing atoms in fixed positions is described with the 
help of multiple wave scattering (MWS) theory. According to it the neutron wave function 
scattered from a single atom at point r = Vi is 

i/)(r) = exp(ikr) — expUkrA- — - — -exp(ik\r — rA), (21) 

\r-ri\ 

where exp(ikr) is the plane wave of the free incident neutron, and b is scattering amplitude, 
which, contrary to the common usage, we define with negative sign. 
If we have N atoms, the neutron wave function becomes 

N b n 

ip(r) = exp(ikr) - ^ Vn ]— 3 — r exp(ifc|r -r n \), (22) 

n=l r r n\ 

where b n is the scattering amplitude of the atom at point r n , and ip n is the amplitude of the 
wave enlightening the n-th atom. The amplitude ip n is a superposition of all the waves coming 
to n-th atom from all the others: 

ip n = exp(ikr n ) - VV] — 3 - exp(ik\r n - r n ,\). (23) 

This is the main system of equations of MWS theory. If we solve it, we find ipn, substitute 
into (|22|) and find the final waves scattered by our set of atoms. 



3.1 Scattering from a single crystalline plane 

We can solve this system, when atoms are identical (bn = b) and arranged regularly on a single 
infinite plane, which we suppose to coincide with coordinate (x, ?/)-plane. For simplicity we can 
imagine the elementary cell of atoms to be a square with side a. From symmetry considerations 
the solution of the system (J23|) is 

ijj n = C exp(ikr n ), (24) 

where r n = an, n = (n x , n y , 0) and n X)V are integers. Substitution of into with b n = b 
gives C = 1/(1 + bS), where S is the sum 

S= exp(ik r n ) eMikrn) , (25) 

which was calculated in jj], and we do not reproduce it here. 

If we substitute (124)1 into (|22|). we obtain wave function scattered by the crystalline plane 

ip(r) = exp(ikr) — Cb^2 e ^^^ Vn ^ exp(ik\r — r n \). (26) 

n \ r ~ r n\ 

We see that C renormalizes b and in the following we omit it. For interpretation we need to 
transform (|26|) . but for that we need some mathematics. 



3.1.1 Some needed mathematics 



First, we need Fourier expansion of the spherical waves: 

d 2 



% f Ob uW 

rjir) = 7^J -^exp^r + ip z \z\), (27) 



where p z = yk 2 — pjj, z axis is selected along normal to the crystalline plane, and vector py 
has components in (x, y) -plane. 

Second, we need to know the representation for a sum of arbitrary numbers f(n) 
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Y^f( m )= jr f f(x)exp(2mMx)dx. (28) 

m=N\ M=—oojy^ 

After substitution of (J2*T|) into (J2SJ) we obtain 

b d ^ 

ip(r) = exp(ifcr) ^exp(ifcna) J — - exp (ip^ (r — an) + ip z \z\^j . (29) 

7T ^ J) z 

With the help of (|2*Hj) the sum is transformed to the other one 

%b f c dp^i 

ip(r) = exp(zfcr) 51 / d 2 nexp(2irmn) exp(ikna) / — -exp Up\\(r — an) + ip z \z\) , 

(30) 

where m = (m x , m y , 0) and m ti!/ are integers. It is easy to integrate (J3U)) over d 2 n and d 2 p, and 
as a result we obtain 

^(r) = e ikr - J2 exp (ik T r + ik ZT \z\) = Q(z < 0) [Mr) + M^} + &(z > 0)^(r), 

y ">ZT 

(31) 

where N 2 = 1/a 2 is the density of atoms on the plane, r = (m x ,m y ,0)/a is a vector of the 



reciprocal plane lattice, k T = k\\ + r, and k ZT = yk 2 — k 2 .. The representation (J3T|) is well 
interpretable. We see that scattering gives diffracted waves. Some of them, e *" ; xT , ||-*fc*T« ) g Q 
back to the space z < 0, and they constitute reflected wave function ip r (r). The others, 
exp fifex?*!! + ik Z TZ^j , go forward to the space z > and they constitute transmitted wave 
function ipt{ r )- Here we use O-functions which are 1 or 0, when inequality in their argument 
are satisfied or not, respectively. 

The set of diffracted waves is infinite, but numerable. The functions ip r {r), ipt{ r ) can be 
represented as infinitely dimensional vectors ip r (f) — exp(ik r r)^f r , ipt(f) — exp(ikr)^ t , where 
ty r j are infinite columns of complex numbers, which can be expanded over orthogonal basis 
|r): ty r j = J2t at^ti-l 7 "). Here a r> t T are expansion coefficients or coordinates of the vectors 
in the basis |r). A basis vector |r) is an infinite column, which contains unity at a position 
corresponding to m x , m y of the vector r, and zero in all other positions. Numeration of 
positions in columns \& can be arbitrary, but the position with m x = m y = we numerate by 
0. The operator k r r has eigen vectors |t) and eigen values k T r\\ — k ZT z. The operator kr has 
the same eigen vectors |r) and eigen values k T r\\ + k ZT z. 

With these notations the state of the incident particle can be represented as exp(fcr)^O) 
where \l/ — \ T — 0), the reflected state is ^ r = r s \l/ , and the transmitted state is = t s ^o, 



where r s and i s = I + f s are reflection and transmission matrices of a single crystalline plane 
with matrix elements 
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Here I is an infinitely dimensional unit matrix, and 5 T y is the Kronecker symbol equal to unity 
for r = r', and otherwise. The total wave function (|3ip can now be represented as 



ip{r) = ®{z < 0) exp(ifcr) + exp(ik r r)r s \l/o + ®(z > 0) exp(ikr)i s ^Q. 



(33) 



We see that diffraction by a single crystal can be reduced to reflection and transmission of a 
one dimensional potential for a particle with an infinitely dimensional spin. 



3.2 Sketch of the dynamical diffraction theory 

Let us imagine the single crystal as a one dimensional system of crystalline planes parallel to 
the entrance surface and separated by distance a. We can consider it as a one dimensional 
potential with symmetric period, consisting of a crystalline plane and two empty spaces of 
width a/2 on both sides of it. Reflection f and transmission t amplitudes of so defined one 
period are easily obtained from r s and t s of a single plane: 

f = Ef a E, i = Ei s E = E 2 + f, (34) 

where E is a diagonal matrix with matrix elements (t\E\t) = exp(ik ZT a/2). This matrix 
describes propagation of the state \I> in empty space. 

First, we can find diffraction from a semiinfinite single crystal. We denote by X the state 
incident upon the second period of the crystal, and put down equations for X and reflection 
matrix R: 

X = i + fRX, R = f + iRX, (35) 

which is similar to ([13)1 . We can resolve the first equation with respect to X and substitute in 
the second one. As result, we obtain 

X = (I-fR)- 1 i R = f + tR(I-fR)- 1 i (36) 

however, solution of this system is more complicated than of (|13|) . because in general matrices 
f, R, X and t do not commute. 

Suppose we could find the solution. Then we are able to find diffraction on a single crystal 
of finite thickness. We define reflection and transmission amplitudes of the crystal 
with iV periods, and put down the system of equations similar to (fT§|): X N = T N + R N RX N , 
R = Rn + TnRX n . We can resolve it with respect to Rn and Tjv, then we obtain 

f N = (I- R 2 )X N [I - RX N RX N }-\ R N = [R- X N RX N ][I - RX N RX N ]~\ (37) 

In the case of Laue diffraction all the intensity goes through the crystal. It means that reflection 
can be neglected, therefore T/v ~ X N ttt N . In the case of Bragg diffraction all the intensity is 
reflected back to the same space, from where the incident wave arrived, therefore X N can be 
neglected and Rn ~ R. 

The simplest case is Laue diffraction. To find the wave function after the crystal we need 
to find eigen function ty q and eigen values X q = exp(zga) of the matrix t: t^ g = \ q ^ q . 



3.3 Some mathematics of matrices 

Our matrix f has a dyadic structure, which means that matrix elements of f have the form 



(Ti\r\Tj) = afij, where a* = -ie^, = 27rbN 2 /k zi , k zi = ^k 2 - (fey + Ti) 2 , = exp(ifc zi a/2), 

6j = e,-. Matrix elements of t are (r.j|t|rj) = e 2 5jj + a^fe,-. 

Let's find eigen vectors and eigen values of t. Eigen vector is J2 x j\ T j)- The equation is 

iJ2xi\Ti) = A^^ilri) — > e 2 Xi + cii^bjXj = Xx^ (38) 

j 

It follows that Xi = diC/(\ — e 2 ), where C = Ej bjXj = C J2j ^j a jl — e f)> from which we get 
the equation for determination of A: 

-iKje 2 _ -2%iN 2 bexp(ik ZT a) _ ^ 
^A-ef ^ k ZT {\ - exp(ik ZT a)) 

If n exponents exp(ik ZTj a) are nearly equal, eq. (J39|) has n solutions Xj close to each other 
and to exp(ik ZT .a). They correspond to n orthogonal eigen states = J2i£n x ji\ T i) > where 
summation is limited mainly to n those Tj, for which exp(ik ZTj a) are nearly the same. 

To find transmitted waves we need to expand the incident wave \l/o over tyj. As a result we 
obtain \l>o = J2jen c j^j- 

Transmission of every eigen vector tyj is equal to exp(iqjNa). Therefore transmitted wave 
function is 

= ^exp^qjNajcj^fj = ^ c, exp^-iVa):^-^), 
from which we see that there are n transmitted waves |tj) with intensity of the transmitted 
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. We obtained the generalization of pandellosung solution for n 



waves \Tj). In the case of n = 2 we have k z0 = k ZT , when \k\\ \ ~ \kn + t\. This is the standard 
case of Laue diffraction. If n — 1, i.e. all exp(ik ZT a) are different, we have refraction with 
q = Jk 2 — -u , where u = 47riV 6, N = N 2 /a = 1/a 3 . 



For Bragg diffraction we need to know R. Let's look at the eq. R = f + tR(I — rR)~ 1 t 
in (}3l)j) . Because of the dyadic form of f, we can find that (1 — fR)" 1 = 1 + CfR, where 
C = 1/(1 — c), and c = J2ij{biRijCLj)- With this in mind we can rewrite eq. for R in the form 

p F 2 rf 2 - ( 1 + E 2 R)f(l + RE 2 ) (oj + e 2 E fc RjkOkHPj + En b n R nj e 2 ) 

1 -Ya.IhA, ' 13 (l-e?e|)(l-E^A) ' 

Matrix elements have a structure of the type Rij = aiXiXjbj/ '(1 — e 2 e 2 ) with amplitudes x^, 
satisfying equations 

1 +Xie 2 J2kXka k bk/(l - e 2 k e 2 ) 
Xi = j (41) 

' ~ Ejfc ajbjXjQkbkXk/ (1 - e\e 2 ) 



We see that the most important are those terms XiXj, for which efe 2 = 1. Let's look at two 
particular cases: specular i — j — 0, and nonspecular j ^ i = Bragg reflections. 

In the specular case, i = j = 0, we denote y = x\, r = ao&o, solve algebraic equation 

V (l " yr 2 /(l - et)) = (l - egyr/(l - e*)) 2 , (42) 
and find R 00 = r y/(l — e^) identical with (fTU|) . where t = + r. 



In the nonspecular case we multiply two eq-s flU)) f° r x o an d Xj, denote y = x Xj, = aibj, 
and obtain the equation 

y (l - yr 0j r j0 /(l - efo 2 )) = (l - e^/(l - e^ 2 )) (l - e>r 00 /(l - efo 2 )) . (43) 

After solution of this equation and substitution Rj = yrj /(l — e^e 2 ) we find for the reflection 
amplitude from state |r = 0) to state |t.,-) the result: 
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ho y // ( 1 + V^Wo) 2 - tootjj -yj(l- y^oTio) 2 - tgotjj ^ 

\ k *Ti \fiX^W^Jof^^j + ^/rojr^f-tootjj 1 

where we substituted ta = ef + r„, = —inieiej, and Kj = 27iN 2 b/k Z T j . The probability of 
such a reflection is 



Wjo = -r- L \Rjo 

&z0 



2 



(45) 



V / ( 1 + V /f D^o) 2_t OO^ + VrbTjo) 2 ~ tootjj 

where the factor k ZT Jk Z Q gives ratio of currents toward the crystal and away from it. 

Because of volume limitations we were not able to present the new dynamical diffraction 
theory in full details. However we hope that from our sketch readers are able to understand 
it. For an exercise we propose to them to answer the question: what will happen if Bragg and 
Laue conditions are satisfied simultaneously? 
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